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ABSTRACT 

Passing stars can perturb the Oort Cloud, triggering comet showers and potentially extinction events 
on Earth. We combine velocity measurements for the recently discovered, nearby, low-mass binary 
system WISE J072003.20-084651.2 (“Scholz’s star”) to calculate its past trajectory. Integrating the 
Galactic orbits of this ~0.15 M 0 binary system and the Sun, we find that the binary passed within only 
52+^ kAU (0.25 +q Q 7 parsec) of the Sun 70 + 4 q kya (ler uncertainties), i.e. within the outer Oort Cloud. 
This is the closest known encounter of a star to our solar system with a well-constrained distance and 
velocity. Previous work suggests that flybys within 0.25 pc occur infrequently (~0.1 Myr -1 ). We 
show that given the low mass and high velocity of the binary system, the encounter was dynamically 
weak. Using the best available astrometry, our simulations suggest that the probability that the star 
penetrated the outer Oort Cloud is ^98%, but the probability of penetrating the dynamically active 
inner Oort Cloud (<20 kAU) is ~10 -4 . While the flyby of this system likely caused negligible impact 
on the flux of long-period comets, the recent discovery of this binary highlights that dynamically 
important Oort Cloud perturbers may be lurking among nearby stars. 

Subject headings: Oort Cloud — Galaxy: kinematics and dynamics — Stars: individual (WISE 
J072003.20-084651.2, Scholz’s star, HIP 85605) 


1. INTRODUCTION 


Perturbations by passing stars on Oort cloud comets 
have previously been proposed as the source of long- 
period c omets visiti n g the planetary region o f the solar 


system (|Oort [1950 


Biermann et al. 1983 Weissman 


1996 Rickman|2014|), and possibly tor generating Earth¬ 


crossing comets th at produce biological extinction events 
(Davis et al.| 1984). Approximately 30% of craters with 
diameters >10 km on the Earth and Moon ar e likely due 


to long-period comets from the Oort Cloud (Weissman 
19961. Periodic increases in the flux of Oort - cloud 


comets due to a hyp otheti cal substc llar companion have 
been proposed (Whitmire & Jackson 1984), however 
recent time series analysis of terrestrial im pact craters 


a re i nconsistent with periodic variations (Bailer-Jones 


2011 ), and sensitive infrared sky surveys have yielded no 


evi dence for an y wide-separation substellar companion 


Hipparcos astrometric data (Perryman et al. 

1997) by 

Garcfa-Sanchez et al. ( 

1999) identified only 

a single 


710; 1.4 Myr in the future at ~0.34 pc), however it 
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is predicted tha t ~12 stars pass within 1 pc of the 
Sun every Myr (Garcfa-Sanche z et al. 2001). A recent 
analysis by Bailer-Joncs (20T3) of the orbits of ^50 ,000 
stars using the revised Hipparcos astrometry from 


van 


Leeuwenj ( |200 7), identified four Hipparcos stars whose 


future flybys may bring them within 0.5 pc of the Sun 
(however the closest candidate HIP 85605 has large 
astrometric uncertainties; see discussion in §3). 

A low-mass star in the solar vicinity in Mono- 
ceros, WISE J072003.20-084651.2 (hereafter W0720 or 
“Scholz’s star”) was recently discovered with a photomet¬ 
ric dista nce of ^ 7 pc and initial spectral classification of 
M9 ± 1 (Scholz 2014). This nearby star likely remained 
undiscovered tor so long due to its combination of prox¬ 
imity to the Galactic plane (b = +2°.3), optical dimness 
(V = 18.3), and low proper motion (~0”.l yr -1 ). The 
combination of proximity and low tangential velocity for 
W0720 (Vtan — 3 kms -1 ) initially drew our attention 
to this system. If most of the star’s motion was radial, 
it was possible that the star m ay have a past or future 
clos e pass to the Sun. Indeed, Burgasser et al.| (|2014 1 
and Ivanov et al. (2014 ) have recently reporte d a high 
positive radial velocity. Burgasser et al. (2014) resolved 
W0720 as a M9.5+T5 binary and provided a trigonomet¬ 
ric parllax distance of 6.0+o '9 pc. Here we investigate the 
trajectory of the W0720 system with respect to the solar 
system, and demonstrate that the star recently (~70,000 
years ago) passed through the Oort Cloud. 


2. DATA & ANALYSIS 

We obtained medium resolution spectroscopy of 
W0720 on UT 17 and 19 November 2013 with the South¬ 
ern African Large Telescope (SALT) and Magellan tele¬ 
scopes, in the o ptical an d nea r -infra red respectively. As 
summarized in Ivanov et al. (2014), we estimated the 
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spectral type of W0720 to be L0±1, and measured a 
radial velocity of 77.6 ±2.5 kms -1 . Recent adaptiv e op¬ 
tics imaging and h igh resolution spectroscopy by |Bur-| 
gasser et al. (2014) indicated that the star is actually a 
low-mass binary classified as M9.5±T5, with projected 
separation 0.8 AU, and a multi-epoch mean radial veloc¬ 
ity of 83.1 ± 0.4 kms -1 . Given the radial velocity accu¬ 
racy, number of reported epochs (11), and imp roved as- 
trometric solution provided by the recent jBurgasser et al. 
(2014) study, we simply adopt their parameters for our 
calculations: a = 07:20:03.21 (ICRS; epoch 2014.0), 5 
= -08°46 , 51".83 (ICRS; epoch 2014.0), )i n = -40.3 ±0.2 
-1 , fis = -114.8 ± 0.4 rnasyr -1 , w = 166 ± 28 mas, 


masyr 

d, 6.0_q.9 PG ^ radial 83.1 ±0.4 km S , Vtotal 

83.2 ±0.4 kms -1 (all la uncertainties). 

Using the Galacti c veloci t y ell ipsoids_and population 
normalizations from |Bensby et al. (2003), we estimate 
that W0720 has probabilities of 77.3%, 22.5%, and 
0 .2% of belonging to the thin disk, thick disk, and 
halo Galactic populations, respectively. The spectral 
classification indices are consistent with a thin disk star 


of ai 
2014 

iproximately solar composition ( 
(classified “dM” based on C = 

Burgasser et al. 
- 1.034 ± 0.018 : 

Lepme et ai. 

2007 

). Classification as a metal-poor 


metallicity and kinematic properties for M-type stars 


( jSavcheva et al.|2014[ ), it is unlikely that W0720 belongs 
to either the thick disk or halo. Among 890 nearby 
solar-typ e stars from a chrom ospheric activity-velocity 
catalog (Jenkins et al. 2011), only 3 have Galactic 


velocity components within ±10 kms -1 of W0720’s 
velocity (HIP 51500, 63851, 117499), and all three are 
less chromospherically active than our Sun (logR’^x 
< -4.95) and have in ferred c hro mos ph eric a ges in the 
range 4-8 Gyr (Mamajek & Hillenbrand 2008). W0720’s 
velocity is similar to” that of the Hercules dynamical 
stream, a kinematic group of stars of heterogenous ages 
and composition, likely perturbed to the solar circle 
from smaller Galactocentric radii d ue to dy n amica l 
interactions with the Galactic bar (Bensby et al. 2007). 
Based on its status as an old thin disk star, consideration 
of the ages of Sun-like stars of similar velocity, and 
the isochronal ages for other Hercules stream members 
( Bensby et al.|2014 ), we adopt an age for W0720 of 3-10 
Gyr (2 ct range). Interpolating the mass-age est i mates 
for the W0720 components from Burgasser et al. (20141 


for solar composition, this age range maps to compo¬ 
nent masses of Mx = 86 ±2 M j up and Ms = 65 ±12 
M ju P (2cr uncertainties). The hyd rogen -burning mass 
limit for stars is near ^7 5 M j up (Saumon & Marley 
2008| Dieterich et al. 2014), hence W0720A is probably 
a low-mass star, and W0720B is probably a brown dwarf. 

We integrated the orbit of W0720 and the Sun 
with a realistic Galactic gravitational potential using 


the NEMO Stellar D ynamics Toolbox (Teuben 1995 
Barenfeld et al. 2013). With the current velocity data, 


we simulated 10 4 orbits of W0720 and the Sun, sam¬ 
pling the observed astrometric values and observational 
uncertainties for W0720 using Gaussian deviates. From 
these 10 4 simulations which take into account the 
observational uncertainties, we find that W0720 passed 
as close as A = 0.252!°;^ pc (lcr; + 0.110 pc 2a) or 


52.0 tlli kAU ( lcr ; ± 22:7 kAU 2cr ) of the Sun. The 
time of closest approach was 70.4+g 4 g 7 (la; j^'g 2a) 
k ±0 Figure 1 shows the distributions of the closest 
approach separations and times for W0720 and the Sun 
for the 10 4 orbit simulations. The median nearest pass 
position from the simulations was (X. Y. Z) = (-0.122, 
±0.120, ±0.185) pc (here quoted in a comoving frame 
centered on the solar system barycenter, where X is in 
the current direction of the Galactic center, Y is in the 
current direction of Galactic rotation, and Z is towards 
the current north Galactic pole), with approximate la 
uncertainties of (ax, ay, az) = (0.047, 0.045, 0.068) 
pc, corresponding to celestial position ( oticRS , Sicrs 
~ 170° ±29°, ±68° ±14°; la unc.), in the vicinity of 
Ursa Major. As a check, we also calculated a linear 
trajectory that ignores the Galactic potential, finding 
that W0720’s closest approach was 70.7 kya at 0.25 
parsec at a barycentric position in Galactic coordinates 
of ( X,Y : Z = -0.125, ±0.119, ±0.185 pc). Note that 
given the short timespan, the predictions of the time 
and position of the W0720-Sun minimum pass from the 
linear trajectory agrees with the more accurate orbit 
integration to better than 2.5%. At its closest pass 
W0720 would have had proper motion exceeding any 
known star: 70” yr -1 (i.e. capable of traversing a full 
moon in 26 years; cf. the current highest proper motion 
star, Barnard’s star, with 10.3"yr -1 ). 


The predicted Galactic coordinates of the nearest 
pass of the binary system (£, b = 135° ± 15°, 47° ± 13°) 
is near one of the two strong peaks in the longitude 
of aphelia distribution of new cl ass I comets from th e 
Oort Cloud (£ = 135° ±15°) (|Matese et al.| |1999|). 
However this appears to be coincidence, as any comets 
on eccentricity ~ 1 orbits from the vicinity of the 
nearest pass ~70 kya would have periods of ~4 Myr, 
and hence require ~2 Myr to reach the inner solar 
system. Also, the two primary peaks in the distribution 
of longitudes of aphelia for long-period comets are 
reasonably explained b y the e ffects of th e sol ar motion 
and Galactic tide (iFeng & Bailer- Jones||2014|). Among 


10 4 simulations, the nearest past separation between 
W0720 and the Sun was A = 0.087 pc (18.0 kAU), 
and this was the only simulation that brought W0720 
within the classical boundary of t he dynamically active 
inner Oort Cloud (“Hills Cloud”) (Hills||1981 |Weissman 
19961 at a < 20 kAU. Approximately 79% of the 


simulations brought W0720 within 0.337 pc of the Sun 
the previously closest estimated pass of a known star 
to the sola r system (G l iese 710; 0.337 pc, 1.4 Myr in 
the future; Garcfa-Sanchez et al. 2001). Approximately 
99.96% of the simulated trajectories brought th e sta r 
within the Sun’s tidal radius of 1.35 pc (Mamajek et al. 
2013), and 98% of the simulated trajectories brought 
the star within the maximal range of semi-major axes 
for retrograde o r biting Oort Cloud comets (^120 kAU) 
(Garcfa-Sanchez et aljl999). The rarity of stellar passes 


with such a small imp act parameter can be a ssesse d 
from the analysis of Garcfa-Sanchez et al. (2001). 
Extrapolating the power-law distribution of minimum 
separations versus cumulative number of encounters per 

9 kya = thousand years ago. 
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Myr from Garcfa-Sanchez et al. (2001), one estimates 
that encounters by stellar systems within 0.25 pc of the 


Sun occur with a frequency of 0.11 Myr -1 or once every 
^9.2 Myr. For comparison, flybys this close (0.25 pc) 
are statistically rare (~2.4%) among encounters by all 
stellar syst ems that penetrate th e Sun’s tidal radius of 
~1.35 pc dMamajek et al7||2013|) , of which ~4.5 occur 
per Myr (Garcia-Sanchez ct ai. 2001). 



Distance of Closest Approach (pc) Number 


Figure 1 . Lower left: Density plot of the distribution of mini¬ 
mum separations between W0720 and the Sun (in pc) versus time 
of minimum separation (in thousands of years ago; “kya”) for 10 4 
Galactic orbit simulations. The maximum semi-major axis for ret¬ 
rograde orbiting Oort Cloud comets is ~0.58 pc (right blue solid 
line), and the dynamically active inner Oort Cloud is at <0.10 pc 
(left blue solid line). The entire region plotted is within the Sun’s 
tidal radius (~1.35 pc). Upper left: Histogram of simulated closest 
W0720-Sun separations (pc). Lower right: Histogram of simulated 
times of closest W0720-Sun separation (kya). 


The flyby’s gravitational interaction is in the regime 
where the Sun’s influen ce on the star’s traj ectory (and 
vice versa) is negligible (Collins & Sari 2010): GM^/bvl 
~ 10 ” 5 ' 6 << 1, where G is the Newtonian constant, 
Mg is the Sun’s mass, b is impact parameter (minimum 
Sun-star separation), and i>* is the velocity of the star 
with respect to the solar system barycenter. Major 
comet showers, where the flux of long-period comets 
increases by factors of >10 are likely limited to cases of 
high mass interlopers passing within ~10 kAU of the 
Sun (Heisler & Tremainejl986), and are exceedingly rare 
(<10 ~ 3 Myr -1 ; extra polating from results by 


Garcfa- 


Sanchez et al. 2001). A perturbed comet with new 


aphelion similar to that of the Sun-W0720 minimum 
separation, but with perihelion in the planetary region 
of the solar system, will have a ~ 26 kAU and P ~ 4.2 
Myr. Although comets throughout the Oort cloud may 
be perturbed by an interloper star onto trajector ies t hat 


would bring them to the inner solar system (Weissman 


19961, the highest density of outer Oort cloud “ comet s 
should exist at smaller orbital distances (Heisler 1990), 
so most of the perturbed comets should originatefrom 
the vicinity of the star’s closest pass to the Sun. Hence, 
any enhancement on the long-period comet flux should 


become manifest ~2 Myr in the future. A proxy indica¬ 
tor of the encounter-induced flux of Oort Cloud comets 


is defined as 7 = M*/v m b (Feng & Bailer-Jones 2014 
using the variables as definecTearlier, and with 1 W* bem; 
the mass of the interloper binary W0720). Using the 10 4 
simulations, we estimate that the W0720 flyby induced 
7 ~ iq - 7 - 48 ± 0.15(1*7) k m ~ 1 sAU“ 1 . Simulations by 


Feng A' Bailer-Jones (2014) suggest that encounters with 
7 < 10 -a 3 are unlikely to generate an enhancement in 
the distribution of longitudes for long-period comets 
compared to that predicted to be generated by Galactic 
tidal effects. All of the 10 4 simulated orbits had 7 < 
10 -, °, hence the pass of the W0720 system should have 
a negligible statistical impact on the flux of long-period 
comets during the coming millenia. 


3. THE CLOSEST KNOWN FLYBY? 

The proximity of W0720’s flyby to the solar system can 
be co mpared to those of H ipparcos stars recently stud¬ 
ied by Bailer-Jones (2014). He identified one star (HIP 
85605) with a flyby closer than the median minimum sep¬ 
aration that we estimated for W0720 (0.25 pc), however 
two future flybys (HIP 89825 [G1 710] and HIP 63721) 
both have A ~ 0.27 pc (0.27 !q; 4 7 pc and 0.27)^'^ pc, 
respectively; 90%CL). |Baiie nJones| (2014) predicts that 
HIP 85605 will pass within A = O.lO+g^g pc (90% CL) 
of the Sun Ml32 k yr in the future. Unfortunately, as 
Bailer-Jones (20141 points out, HIP 85605 is a visual bi- 
nary, and there is considerable dispersion in its published 
proper motions and its Hipparcos parallax is of low ac¬ 
curacy. The solution which brings HIP 85605 close to the 
Sun relies on the Tycho-2 pr oper moti o n ( p a , ps = +4.0, 


-7.6 ± 2.0, 1.9 rnasyr 4 ; Hpg et al. 2000), the revised 


Hipparcos parallax of w = 146.84 + 29.81 mas, and the 


Pulkovo radial velocity of -21.0 + 0.3 km/s (Gontcharov 
20061 as reported in the XHIP compiled catalog; |Ander- 
son"& Francis|(2012)). If any of these values are substan- 


( 2012 )). it any 
tiaily in error, then the proximity of the flyby solution 
is likely to be spurious. Both t he original and revised 


Hipparcos astrome tric solutions (Perryman et al. 1997 
van Leeuwen|2007 ) ha d much larger prope r motions than 


the Tycho -2 solution. Bailer-Jones (20141 mentions that 


a reanalysis of the Hipparcos astrometric solution by F. 
van Leeuwen suggests u that the Hipparcos-^ parallax of 
Hip 85605 and its (relatively large) uncertainty are valid, 
but that on account of the large residuals and the complex 
nature of this system the solution should be treated with 
caution .” 

Unfortunately, the Hipparcos parallax for HIP 85605 


leads to som e astrophysical inconsistency. As Bailer- 
Jones (2014) note, HIP 85605’s B-V color (1.1 mag) is 


consistent with a K-type star, whereas the absolute mag¬ 
nitude calculated using its revised Hipparcos parallax 
(My = 11.86 + 0.45 mag) is more consistent with an M 
dwarf. We find that HIP 85605’s colors from Hipparcos 
and 2MASS (B-V = 1.10 ± 0.11, V-K s ~ 2.66 ± 0.07, J-H 
~ 0.53 ± 0.03, H-K s _c; 0.12 ± 0.03)are all consistent with 
a ~K4 dwarf (Pecaut & Mamajek] 2013| , in agreement 
with the spectr al template fitting analysis of Pic kles fe| 
Depagne (2010, ; which also yielded K4V). Dave Latham 
(priv. comm.) has visually confirmed that five spectra 
of HIP 85605 taken with the CfA Digital Speedometer 
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(as reported in 

Garcfa-Sanchez et al. 

1999) are indeed 

the National Research Foundation (NRF) of 

South 

consistent with a typical K dwarl. As I 

HP 85605 is spec- 

Africa. We thank Dave Latham, Alice Quillen, 

Kevin 


trophotometrically a K dwarf, the Hipparcos- derived ab¬ 
solute magnitude (My ~ 11.9) places the star nearly 5 
magnitudes fainter than the main sequence (M y S ~ 7.0) 
for a star of its color - too faint to be simply a metal- 
poor dwarf, yet too red and luminous for a white dwarf 13 ) 
Hence, most likely HIP 85605 is a ~K4V at d ~ 60 pc and 
the Hipparcos parallax is erroneous. At d ~ 60 pc, HIP 
85605 has a slightly different velocity (U,V,W = -10.7, 
-14.3, -11.6 kms -1 ) compared to that used in the simu¬ 


lation by Bailer-Jones (2014), and its “flyby” of the Sun 
correspondingly moves further into the future and further 
away (~2.8 Myr in future, ~10 pc away). We conclude 
that HIP 85605 is unlikely to penetrate the Sun’s Oort 
Cloud, and that W0720 now appears to have the closest 
flyby of any known star. 

4. DISCUSSION 


Given its current visual magnitude of V ~ 18.3 (Ivanov 
et al.|2014|), at its closest approach of ~0.25 pc, Schoiz’s 
star ( W0720) would have had an apparent magnitude of 
V — 10.3, brighter than the current nearest star (Prox- 
ima; V = 11.2) but still much dimmer than the faintest 
naked eye stars ( V ~ 6). However, W0 720 is an ac¬ 
tive M dwarf star (Burgasser et al. 2014), and E-band 
flares have been observed to exceed 9 mag on timescales 
of minutes amongst such s tars, and brief flares e xceeding 
12 mag may be possible (Schmidt et al. 20141. Flares 
amongst the coolest M dwarfs have been witn essed with 
energies of ~10 34 ergs (Schmidt et al. 2014) and lumi¬ 
nosities of ~10 29 -10 3 ° ergs -1 . If W0720 experienced 
occasional flares similar to t hose of the active M8 star 
SDSS J022116.84+194020.4 flSchmidt et al.l|2014|), then 
the star may have been rarely visible with the naked eye 
from Earth (V < 6; AE < —4) for minutes or hours dur¬ 
ing the flare events. Hence, while the binary system was 
too dim to see with the naked eye in its quiescent state 
during its flyby of the solar system ~70 kya, flares by 
the M9.5 primary may have provided visible short-lived 
transients visble to our ancestors. 

Improved astrometry f or W0720 via groun d-based 


telescopes (Dieterich et al. 2014) 


via 
or Gaia 


(de Bruijne 


20121 , and further radial velocity monitoring, should help 


reduce the uncertainties in the flyby timing and mini¬ 
mum separation between Schoiz’s star and the solar sys¬ 
tem. Past systematic searches for stars with close flybys 
to the solar system have been understand ably focused 
on the Hipparcos astrometric catalog (Garcfa-Sanchez 
et al.||1999| Bailer- Jones||2014 ), however it contains rela¬ 
tively few M dwarfs relative to their cosmic abundance. 
Searches in the Gaia astrometric catalog for nearby M 
dwarfs with small proper motions and large parallaxes 
(i.e. with small tangential velocities) will likely yield ad¬ 
dition candidates. 
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